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Preface 


This  thesis  represents  a small  part  of  the  Air  Force 
Flight  Dynamics  Lab's  effort  to  improve  the  stopping 
performance  of  aircraft.  Major  impetus  was  provided  by 
the  sensitivity  study  of  the  stopping  parameters  done  by 
the  Boeing  Commercial  Aircraft  Company  under  the  auspices 
of  the  Combat  Traction  program  in  1974.  A major  finding 
of  this  program  was  the  importance  of  tire  performance  in 
the  stopping  equation.  This  thesis  is  intended  to  improve 
tire  performance  by  improving  the  procedure  under  which 
tires  are  designed  and  tested. 

In  this  study  I have  employed  the  analog  computer  to 
solve  the  aircraft  equations  of  motion  and  to  predict  main 
gear  loads  of  an  aircraft  under  varying  conditions  of  gross 
weight,  touchdown  velocity,  sink  rate  and  runway  condition. 
To  do  this  I have  comgined  the  efforts  of  Boeing  in  their 
Brake  Control  Simulator,  McDonnell  in  their  Motion  Base 
Simulator,  and  NASA  Langley  Research  Center  in  their  study 
of  tire  traction.  The  result  is  a three  degree  of  freedom, 
pure  analog  system  simulation. 

Major  thanks  go  to  Jerry  Schumacher  of  AFFDL,  Mecha- 
nical Branch  for  providing  valuable  insight,  contacts,  and 
information  essential  to  the  successful  completion  of  this 
effort.  Major  thanks  also  go  to  my  thesis  advisor,  Dr. 
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Robert  A.  Calico  for  his  guidance  and  assistance  with  the 
analog  computer,  and  to  my  thesis  committee  for  their 
constructive  criticism.  Finally,  I would  like  to  thank 
Barbara  Barnes  for  her  assistance  in  editing  and  typing 
and  my  wife,  Barbara,  for  her  patience,  understanding, 
and  support. 


John  A.  Skorupa 
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Abstract 


Main  gear  load  vs.  time  is  predicted  for  an  F-4E  air- 
craft through  the  use  of  subsystem  modeling  and  analog 
computation.  Subsystems  modeled  are  the  aerodynamics, 
engine  dynamics,  vertical  strut  dynamics,  fore-aft  strut 
dynamics,  tire/wheel  dynamics , brake  dynamics,  and  anti- 
skid dynamics.  The  problem  is  restricted  to  a landing 
sequence  with  three  degrees  of  freedom  permitted  for  the 
aircraft.  Aerodynamics  are  based  on  constant  coefficients 
of  lift,  drag,  and  pitching  moment.  A drag  chute  is  also 
employed.  Engine  dynamics  are  based  on  a linear  thrust  vs. 
velocity  schedule.  The  strut  dynamics  are  modeled  by  a 
mass-spring-damper  system.  The  tire/wheel  dynamics  sub- 
system applies  Newton's  Second  Law  to  derive  the  wheel  ve- 
locity, and  calculates  the  wheel  slip  ratio  and  ground-tire 
coefficient  of  friction.  Brake  dynamics  are  based  on  a 
schedule  of  brake  torque  vs.  brake  pressure.  Antiskid  dyna- 
mics model  the  Hytrol  Mark  II  antiskid  system.  Stopping 
distances  from  simulation  are  compared  to  flight  test  data 
to  verify  the  model. 

Results  from  the  simulation  agree  with  flight  test 
data.  A schedule  of  main  gear  load  vs.  velocity  is  pro- 
posed as  an  alternative  to  current  tire  testing  practice. 
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SYSTEM  SIMULATION  IN  AIRCRAFT 
LANDING  GEAR  AND  TIRE  DEVELOPMENT 

I.  Introduction 


In  a theoretical-experimental  vein,  this  thesis  will 
explore  one  aspect  of  landing  gear  and  tire  developmental 
practice  with  the  purpose  of  proposing  an  alternative 
procedure  based  on  experimental  findings.  The  approach 
will  be  to  divide  the  aircraft  into  a group  of  subsystems 
and  determine  the  net  effect  of  their  interaction. 

Background 

Landing  gear  and  tire  development  is  critically 
concerned  with  the  operating  environment  of  the  aircraft. 
For  military  aircraft  this  environment  is  determined  by 
a series  of  military  specifications  which  spell  out  per- 
formance goals.  These  mil  specs,  as  they  are  commonly 
referred  to,  supposedly  take  into  consideration  the 
operating  limits  of  the  various  aircraft  systems  which 
interface  with  the  landing  gear  and  tires.  Historically, 
mil  specs  for  the  design  of  tires  have  set  performance 
standards  based  on  static  analysis  of  forces  projected 
through  a velocity  vs.  distance  schedule  as  shown  in  Fig. 
1.  The  glaring  omission  with  these  is  the  absence 
of  dynamic  effects.  For  example,  consider  a brake  system 
dependent  on  antiskid  system  inputs.  If  the  two  were 
designed  without  considering  brake  pressure  lag  or  the 
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Load/Wheel 
(xlO^  lbs)  18 


Velocity 
(xlO  Knots) 


Velocity 


Runway  Distance  (xlCr  ft.) 


Fig.  1.  Military  Specification  LoadAelocity  Schedule 
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frequency  response  of  the  antiskid  system,  it  might  be 
possible  to  have  the  two  systems  operating  out  of  phase 
such  that  braking  never  occured!  Applying  this  analogy  to 
landing  gear,  the  test  pilot  of  a prototype  aircraft 
developed  under  existing  mil  specs  for  tires  and  landing 
gear  has  little  assurance  that  dynamic  effects  will  not 
interact  catastrophically  on  the  first  landing  or  aborted 
take  off. 

Problem  Statement 

A more  logical  approach  to  landing  gear  design  would 
be  to  model  interacting  aircraft  systems  on  an  analog  com- 
puter so  that  realistic  design  specifications  could  be 
based  on  dynamic  simulation  of  expected  forces.  A schema- 
tic diagram  of  how  this  approach  could  work  is  presented 
in  Pig.  2 . The  analog  computer  simula4Mfc  - ' aero- 
dynamics, engine  thrust,  vertical  strut  dynamics,  tire/ 
wheel  dynamics,  braking  dynamics,  antiskid  dynamics,  and 
fore-aft  strut  dynamics  would  combine  to  provide  either 
real-time  load/velocity  vs.  distance  information  for 
testing  or  a schedule  of  load/velocity  vs.  distance  for 
design  specifications.  The  problem  of  dynamic  force 
simulation,  therefore,  is  achievable  if  the  modeling  of 
the  aircraft  dynamics  and  pertinent  subsystems  shown  in 
Fig.  2 is  performed. 
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Scope  of  Effort 

As  a proof  of  principle  this  analysis  will  be  limited 
to  studying  an  aircraft  in  a landing  sequence,  operating 
in  ground  effect.  To  facilitate  gathering  data,  this 

analysis  will  be  restricted  to  an  aircraft  for  which  well 

j 

accumulated  data  exists,  namely  an  F-4E.  In  particular, 
the  F-4E  will  be  considered  to  have  gear  down  with  no 
external  stores,  leading  and  trailing  edge  flaps  deployed 
and  a boundary  layer  control  system  operating.  Since  the 
landing  dynamics  are  constrained  to  a runway  environment 
with  limited  lateral  travel,  only  three  degrees  of  freedom 

I 

will  be  assumed:  forward  translation,  vertical  transla- 

tion, and  pitching  rotation.  Despite  this  latter  assump- 
tion, landing  gear  compression  limits  for  the  F-4E  also 
warrant  the  assumption  of  constant  angle  of  attack  on  the 
runway.  Pitching  moment  due  to  idle  thrust  will  be  neglec- 
ted. Pilot  inputs  will  also  be  neglected;  the  controls 
will  be  fixed  with  brakes  being  applied  two  seconds  after 
touchdown  and  drag  chute  deployed  four  seconds  after  touch- 
down. The  aircraft  will  be  assumed  to  be  a rigid  body. 

Landing  gear  dynamics  will  be  assumed  to  be  analogous  to 

mass-spring-damper  dynamics  with  the  mass  of  the  wheel  and 

strut  concentrated  in  the  wheel  and  a massless  spring- 

damper  combination  simulating  the  strut.  Small  angle  static 

deviation  of  the  struts  from  vertical  will  be  neglected.  j 

l I ^ 

I Tire  deformation  will  also  be  neglected,  so  that  vertical 

gear  loads  will  act  through  the  center  of  the  wheel. 
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Proposed  Procedure 

The  modeling  of  this  landing  sequence  can  be  broken  in- 
to two  related  sections.  The  first  section  is  concerned 
with  producing  the  main  gear  load  and  aircraft  velocity. 

This  section  would  model  the  aircraft  dynamics  if  an  average 

. 

ground-tire  coefficient  of  friction  is  assumed.  The  second 
section  is  concerned  with  calculating  a variable  ground- 
tire  coefficient  of  friction.  The  aircraft  subsystems  that 
are  primarily  concerned  in  the  generation  of  gear  loads  and 
velocity  are  the  aerodynamics,  the  engine  dynamics,  and 
the  vertical  strut  dynamics. 

The  aerodynamics  subsystem  produces  lift,  drag,  and 

I ' 

pitching  moment,  which  are  each  equal  to  a constant  times 
the  square  of  the  velocity  for  constant  angle  of  attack, 

* altitude,  and  wing  area.  The  velocity  required  for  these 

calculations  is  found  by  solving  the  aircraft  equations  of 
motion.  The  acceleration  is  integrated  to  obtain  the  velocity. 
, The  engine  dynamics  are  also  only  a function  of  velo- 
city since  a constant  altitude  is  assumed  and  idle  thrust 
is  selected  prior  to  touchdown.  At  low  speeds  the  rela- 
tionship is  essentially  linear,  so  that  (thrust)  *=  (static 
idle  thrust)  + (K) (velocity) , where  K is  a constant. 

The  vertical  strut  dynamics  simulation  will  be  such 
that  a change  in  the  gear  load  is  produced.  Since  the 
aircraft  dynamics  is  summing  forces  and  accelerations  in 
each  of  the  three  directions  corresponding  to  the  three 
degrees  of  freedom,  a static  gear  load  is  already  known, 
given  the  geometry  of  the  aircraft.  The  change  in  gear 
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load  will  be  added  to  the  static  gear  load  to  correct  for 
dynamic  effects.  The  inputs  to  the  strut  dynamics  required 
to  calculate  the  change  in  gear  load  are  the  vertical  velo- 
city and  displacement  at  the  strut.  Since  vertical  transla- 
tion and  pitching  rotation  are  permitted,  the  strut  vertical 
velocity  and  displacement  will  both  have  components  due  to 
vertical  movement  and  rotational  movement. 

Calculating  the  variable  ground-tire  coefficient  of 
friction  requires  modeling  the  tire/wheel  dynamics,  the 
braking  dynamics,  the  antiskid  dynamics,  and  the  fore-aft 
strut  dynamics.  The  antiskid  subsystem,  basically, 
senses  wheel  speed  to  determine  a skid  and  to  vent  brake 
pressure  when  one  is  detected.  The  output  is  a modified 
brake  pressure  with  wheel  speed  as  an  input.  The  wheel 
speed  is  determined  by  integrating  the  wheel  acceleration 
produced  by  the  difference  in  brake  torque  and  ground 
torque . 

The  brake  dynamics  subsystem  produces  a brake  torque 
from  a known  schedule  of  brake  torque  vs.  brake  pressure. 
This  torque  is  then  corrected  for  the  dynamic  effects  of 
fade,  pressure  lag,  and  torque  peaking  which  are  functions 
of  brake  pressure,  wheel  speed,  brake  torque,  and  antiskid 
cycling. 

The  tire/wheel  dynamics  subsystem  calculates  a wheel 
slip  ratio  by  comparing  the  wheel  velocity  after  braking 
with  the  velocity  of  the  aircraft.  This  value  is  matched 
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against  a preprogrammed  schedule  to  determine  the  ground- 
tire  coefficient  of  friction.  This  coefficient  of  friction 
is  then  supplied  to  the  aircraft  dynamics  model  to  provide 
a corrected  braking  force  in  the  foreward  translation 
equation. 

The  fore-aft  strut  dynamics,  once  again,  employ  a 
mass-spring-damper  analogy  to  produce  an  oscillatory 
velocity  that  modifies  the  wheel  speed  sensed  at  the  run- 
way interface.  The  inputs  required  are  a coefficient  of 
friction  and  a main  gear  load.  The  coefficient  of  fric- 
tion is  produced  by  the  tire/wheel  dynamics,  and  the  main 
gear  load  is  produced  by  the  vertical  strut  dynamics  and 
the  aircraft  dynamics. 

The  simulation  used  three  Electronics  Associates, 

Inc.,  TR-48  analog  computers.  Landing  distance  produced 
by  the  simulation  will  be  compared  to  flight  test  stopping 
distances  to  verify  the  model.  Load/velocity  vs.  distance 
schedules  will  then  be  prepared  and  compared  to  military 
specifications.  The  net  result  of  this  system  simulation 
will  be  a realistic  schedule  of  load/velocity  vs.  distance 
that  can  serve  as  a military  specification  for  tire 
development  and  testing.  Since  all  of  the  data  required 
to  complete  this  model  is  easily  generated  early  in  the 
design  sequence,  this  method  could  serve  as  a new  procedure 
in  landing  gear  development. 
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II.  Procedure 


Introduction 

The  procedure  used  herein  is  divided  into  four  basic 
areas  of  concern:  (1)  force  analysis  and  subsystem  desig- 

nation; (2)  subsystem  modeling  and  total  system  integra- 
tion; (3)  data  acquisition;  and  (4)  model  verification. 
Assumptions  made  will  be  discussed  followed  by  formulation 
of  mathematical  equations  of  motion  and  conversion  to 
computational  equations. 

Force  Analysis 

To  achieve  the  benefits  of  subsystem  modeling  a natu- 
ral starting  point  is  with  an  analysis  of  the  forces 
acting  on  the  aircraft.  For  an  aircraft  in  a landing  se- 
quence, an  appropriate  choice  of  the  coordinate  system  will 
aid  considerably.  For  example,  selecting  the  x-direction 
to  be  in  the  direction  of  horizontal  motion,  as  opposed  to 
the  longitudinal  axis  of  the  aircraft,  and  the  z-direction 
to  be  perpendicular  to  the  x-direction  and  pointing  toward 
the  center  of  the  earth,  forces  such  as  lift  and  drag  are 
easily  reconciled.  The  remaining  y-direction  corresponds 
to  a right  hand  rule  Cartesian  coordinate  system  as  de- 
picted in  Fig.  3. 

By  recognizing  that  the  no-crosswind  landing  sequence 
properly  performed  involves  no  rotation  about  the  z-axis 
or  x-axis,  nor  any  motion  side-ways,  it  is  reasonable  to 
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assume  that  the  corresponding  moments  and  forces  are 
perpetually  balanced  and  are  of  no  significance  to  this 
problem.  In  specifying  this,  however,  several  relevant 
assumptions  are  implied  and  should  be  stated. 

• There  can  be  no  crosswind  component. 

• No  forces  should  be  transmitted  to  the  tires  with 
components  in  the  y-direction,  nor  should  any 
moments  be  transmitted  to  the  tires  about  the  x- 
axis  or  z-axis. 

• The  yaw  and  bank  angles  are  zero. 

« Symmetrical  runway  conditions  exist  on  either  side 
of  aircraft  centerline. 

With  these  assumptions  a symmetrical  landing/roll-out 
condition  is  specified  and  the  six  equations,  resulting 
from  summing  forces  and  moments  about  the  aforementioned 
axes,  uncouple  and  reduce  to  three  equations. 
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EFy  = 0 Lateral-Directional  Equations 

£MX  = 0 (Perpetually  balanced,  not  rele- 


EMZ  = 0 vant  with  respect  to  assumptions) 


To  further  define  the  forces  it  is  necessary  to  know 
something  about  the  type  of  aircraft  and  its  configuration. 
For  this  study  the  McDonnell  F-4E  was  selected,  based  on 
availability  of  data  and  future  applicability  of  potential 
results.  By  referring  to  the  F-4E  we  specify:  a single 

< nose  gear  strut  with  two  unbraked  tires,  two  main  gear 

struts  with  single  braked  tires  each;  leading  and  trailing 
edge  flaps;  a boundary  layer  control  system;  the  Hytrol 
Mark  II  antiskid  system;  a drag  chute;  and  two  General 
Electric  J-79-17  engines  with  thrust  inclined  from  the 
horizontal  at  an  angle,  4»t  = 5.25°.  By  further  assuming 
that  the  landing  simulation  begins  with  the  tires  in 
contact  with  the  runway,  it  is  now  possible  to  locate 
force  components  on  a profile  of  the  selected  aircraft 
(Fig.  4). 

By  lumping  related  forces  together,  several  distinct 
groupings  form.  The  lift,  drag,  and  pitching  moment  all 
are  aerodynamic  in  origin  and,  thus,  are  considered  to  be 
a product  of  the  aerodynamic  subsystem.  The  engine  thrust 
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Body  Force  Diagram 


is  unique  and  is  considered  to  be  a subsystem  in  itself. 

The  normal  forces  on  the  landing  gear  arise  from  two  sepa- 
rate mechanisms.  The  first  is  the  load  caused  by  the  weight 
of  the  aircraft  as  the  lift  decays;  the  second  is  the  load 
caused  by  touchdown  impact  and  dynamic  force  interaction. 

These,  then,  are  grouped  together  and  referred  to  as  the 
vertical  strut  dynamics.  If  it  is  assumed  that  the  fuel  con- 
sumption at  idle  thrust  and  tire  surface  losses  during  the 
landing  are  so  small  as  to  be  insignificant  when  compared 
to  the  mass  of  the  aircraft,  it  can  then  be  assumed  that  the 
weight  is  a constant.  This  being  so,  no  separate  subsystem 
is  required  for  modeling  purposes.  The  braking  and  rolling 
friction  forces,  likewise,  are  similar  in  that  they  both  are 
dependent  on  the  normal  force  of  their  respective  struts  and 
the  respectively  different  coefficients  of  friction.  For 
the  nose  gear,  the  coefficient  of  rolling  friction  is  assumed 
to  be  a constant.  For  the  main  gear,  the  coefficient  of 
braking  friction  is  dependent  on  the  ratio  of  wheel  slip 
speed  to  ground  speed.  The  computation  of  this  coefficient 
then  is  sufficiently  complicated  to  warrant  a tire/wheel 
dynamics  subsystem.  This  tire/wheel  dynamics  subsystem, 
however,  is  not  independent  since  it  must  rely  on  the  wheel 

j 

speed.  Since  the  wheel  soeed  is  rarely  equal  to  the  aircraft 
speed,  owing  to  braking  action,  a brake  dynamics  subsystem  is 
required.  If  a brake  dynamics  subsystem  is  considered,  the 
antiskid  system  must  also  be  accounted  for.  Considering 
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the  geometry  of  the  main  gear  strut  and  the  importance  of 
wheel  speed  to  the  calculations,  one  recognizes  a component 
of  wheel  speed  resulting  from  fore-aft  strut  oscillation. 

To  account  for  this  component  of  wheel  speed,  a fore-aft 
strut  dynamics  subsystem  is  required.  If  we  assume  the 
surface  of  the  earth  serves  as  an  inertially  fixed  reference 
point,  considering  the  speeds  involved,  then  these  subsystems 
(aerodynamics,  engine,  vertical  strut,  tire/wheel,  brake, 
antiskid,  and  fore-aft  strut)  plus  weight  of  the  aircraft 
account  for  all  the  external  forces  acting  on  the  aircraft. 

If  the  only  requirement  was  to  find  the  total  load  on  the 
nose  and  main  gear  combined,  the  number  of  unknowns  would 
equal  the  number  of  independent  equations;  however,  the 
ultimate  goal  is  to  know  the  main  gear  load,  thus,  requir- 
ing an  additional  equation.  This  equation  can  be  found  by 
considering  the  static  loads  on  the  gear  struts.  Since 
this  proportion  is  0.1396  to  0.8604,  nose  gear  load  to  main 
gear  load  for  a .33c  center  of  gravity  location,  the  problem 
is  now  solvable.  The  actual  solving  of  these  four  equations 
will  represent  the  integration  of  all  the  subsystem  models 
and  be  referred  to  as  the  aircraft  dynamics. 

Subsystem  Modeling 

Aerodynamics.  For  the  landing  problem,  the  altitude 
is  a constant  and,  since  planform  changes  do  not  occur  for 
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the  F-4E  from  touchdown  to  stop,  the  wing  area  and  chord 


length  are  also  a constant.  If  the  strut  compression 
limits  for  the  F-4E  are  considered,  one  observes  that 
from  static  equilibrium  the  main  gear  travel  is  plus 
or  minus  one  inch,  while  the  nose  gear  travel  is  plus 
or  minus  12  inches  {Ref  8:V2,  V5) . Taking  the  worst 
combination  of  the  two  extremes  acting  over  the  wheel 
base  results  in  an  angle  of  attack  change  of  plus  or 
minus  2.66°.  Since  this  will  most  likely  occur,  if 
ever,  only  at  touchdown  and  will  be  damped  out  rapidly 
by  viscous  damping  action  of  the  struts,  its  effect 
over  the  duration  of  a landing  roll  will  be  minimal. 
i Therefore,  it  is  assumed  that  the  angle  of  attack  on 

the  runway  is  a constant,  equal  to  the  angle  of  inci- 
dence of  the  wing.  By  assuming  this,  it  follows  that 
the  coefficients  of  lift,  drag,  and  pitching  moment  are 
a constant  for  a given  configuration. 

In  terms  of  configuration,  a normal  landing  calls 
for  leading  edge  flaps  set  to  0°/60°/60°  running  from 
inboard  to  outboard,  trailing  edge  flaps  set  to  60°,  no 
speed  brake,  no  spoilers,  and  no  thrust  reversing  (Ref 
11:4.2).  Furthermore  to  obtain  the  correct  lift,  drag, 
and  pitching  moment,  the  coefficients  must  be  based  on 
gear  down  configuration  with  a correction  for  being  in 
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ground  effect.  Implied  in  this  assumption  of  constant 
configuration  is  the  assumption  that  the  flight  controls 
are  set  at  a constant  position,  untrimmed.,  with  no  pilot 
inputs.  The  rudder  and  ailerons  are  set  to  neutral  by 
the  assumption  of  no  lateral-directional  forces.  The  ele- 
vator is  set  to  correspond  to  a slight  stick  aft  position 
for  which  data  exists — nominally  minus  eight  degrees.  The 
drag  polar  that  provides  Cl,  Cq,  and  Cm  corresponding  to 
the  above  conditions  (Ref  11:5.4)  sets 


ClBASIC  °*195 


CDBASIC  = °*140 


’BASIC 


0.025 


The  correction  for  ground  effect  based  on  6g  ■ -8°,  and 
angle  of  incidence  of  1°  provides  the  following  equations 
(Ref  11:14.2) 


CL  = CL 


BASIC 


CD  * CD 


BASIC 


+ 0.077  = 0.272 


- 0.023  = 0.117 


CM  " CMBASIC  - 0-004  = 0.021 


The  remaining  term  in  the  aerodynamic  subsystem  is 
the  drag  and  pitching  moment  contribution  of  the  drag 


chute.  For  the  16  ft.  diameter  ring-slot  chute  carried  by 
the  F-4E,  the  drag  coefficient  is  0.1875  (Ref  9:7.12.2). 
This  implies  that 


Ddc  = (0.1875) (1/2) pSV' 


The  pitching  moment  is  derived  by  multiplying  D<jc  by  the 
moment  arm  from  the  center  of  gravity  to  the  vertical 
height  of  the  drag  chute  attachment  point. 

For  computation  of  landing  distances,  the  flight  test 
data  is  based  on  brakes  being  applied  two  seconds  after 
touchdown  and  drag  chute  fully  deployed  four  seconds  after 
touchdown  (Ref  9:7.12.2).  Drag  chute  initiation  is  simul- 
taneous with  brake  application  but  requires  two  seconds 
deployment  time.  Drag  chute  effectiveness  is  assumed  to 
be  a step  input  based  on  pilot  critique  of  the  McDonnell 
Motion  Base  Simulator's  one  second  ramp  fade-in  procedure 
(Ref  3:58). 

The  only  unknown  in  the  four  aerodynamic  equations 
is  V,  the  aircraft  true  airspeed.  This  is  provided  during 
the  process  of  solving  the  equations  of  motion. 

Given  the  physical  values  for  the  terms  in  the  aero- 
dynamic equations  as  listed  in  Appendix  A,  it  is  now  possi- 
ble to  arrange  the  equations  for  the  analog  computer. 
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L = 0.1714x^  lbs. 


D = 0.0737X2  lbs. 


Mp  = 0.2123x2  ft-lbs, 


Ddc  = 0.1182X2 


zdcDdc  = 0.4302x2 


To  scale  these  equations  for  the  analog  computer  it  is 
necessary  to  know  the  maximum  value  that  the  variables  can 
have.  Assuming  x has  a maximum  value  of  280  ft. /sec. , 
the  equations  are  scaled  in  terms  of  x as  a fraction  of 
its  maximum  value  and  take  the  form 


[L/l]  = 0.1714 [x/280]22802  = 13437 . 76 [x/280] 2 


[D/1]  = 0.0737 [x/280]22802  = 5778.08 [x/280]2 


[Mp/1]  = 0.2123 [x/280]22802  = 16644 . 32 [x/280] 


tDdc/U  = 0.1182  [x/280]22802  = 9266 . 88  [x/280] 


[ZdcDdc/H  » 0.4302 [x/280] 22802  = 33,727.68 [x/280]2 
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Engine  Dynamics . Thrust  from  the  two  GE  J-79-17 
engines  at  idle  power  is  not  a constant,  but  varies  with 
velocity.  In  the  context  of  the  landing  sequence,  the 
thrust  varies  linearly  with  velocity  from  touchdown  to 
stop,  assuming  a maximum  landing  velocity  of  M = 0.25 
(Ref  9:5.1.38).  At  velocity  of  M = 0 the  static  idle 
thrust  of  470  lbs. /engine  is  cited.  At  M = 0.2,  the  thrust 
is  only  50  lbs. /engine.  Therefore,  assuming  a linear 
relationship,  single  engine  thrust  takes  the  form 


T 


(50-470) 

(0.2-0)  fa) 


x + 470  lbs 


where  (a)  represents  the  speed  of  sound  at  sea  level  or, 
nominally,  1115  ft/sec.  This  being  the  case,  for  two 
engines  the  thrust  equation  takes  the  form 

T * -3.7668 x + 940  lbs. 

Scaling  this  equation  for  the  analog  computer  requires 
the  recognition  of  the  maximum  value  for  x.  Recalling  that 
280  ft. /sec.  was  used  in  the  aerodynamics  section,  the 
equation  takes  the  form 

[T/1J  = -1054.704 [x/280]  + 940  lbs. 
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Vertical  Strut  Dynamics . For  controlability  reasons 
the  impact  of  a landing  aircraft  with  a runway  must  be 
such  that  the  aircraft  remains  in  contact  with  the  runway 
throughout  the  period  from  touchdown  to  stop.  Accomplishing 
this  requires  that  the  landing  gear  absorb  the  landing 
energy  and  dissipate  it  over  a period  of  time.  Usually 
this  is  done  through  the  combined  use  of  a compressible 
fluid,  acting  as  a spring  force,  and  an  incompressible 
fluid,  acting  as  a viscous  damper.  Each  landing  gear, 
therefore,  is  modeled  as  a mass-spring  damper  system  as 
shown  in  Fig.  5,  with  two  degrees  of  relative  freedom 
allowed  for  the  aircraft. 


Mass  of  Aircraft 


Nose  Gear 
_i_ 


J 


Main  Gear 


Fig.  5.  Vertical  Strut  Analogy 


! 
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In  applying  this  concept  to  the  landing  aircraft,  the 
force  of  the  aircraft  on  the  gear  can  be  assumed  to  have 
two  components,  a quasi  static  load  contribution  and  a 
dynamic  load  contribution.  The  quasi  static  load  contri- 
bution is  the  first  mechanism  described  in  the  force 
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analysis  section  whereby  the  landing  gear  is  loaded  and  it 
accounts  for  all  the  external  forces  on  the  aircraft 
acting  through  the  center  of  gravity.  These  are  distribu- 
ted with  respect  to  the  location  of  the  aircraft  center  of 
gravity  and  the  landing  gear.  The  result  is  86.04%  of  the 
load  is  on  the  main  gear  and  13.96%  of  the  load  is  on  the 
nose  gear.  The  dynamic  load  contribution  is  the  second 
mechanism  and  stems  from  the  vertical  and  rotational  motion/ 
displacement  of  the  aircraft.  With  this  in  mind  the  forces 
on  the  gear  become 


smg  Fsmg  + F<img 

* <0-8604  load)  + + Ksm9zsmg) 


Fzng  Fsng  + Fdng 

" (0*1396  load)  + (Kdng*sng  + Ksngzsng> 


In  these  equations  the  (load)  term  represents  the  sum  of 
all  the  vertical  loads;  lift,  weight  and  engine  thrust. 
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To  maintain  a constant  attitude  on  the  runway  for 
changing  gross  weights,  the  nose  gear  and  main  gear  are 
designed  to  displace  equally  for  changing  external  forces 
acting  through  the  aircraft  center  of  gravity.  Since  the 
static  gear  loads  result  from  forces  acting  through  the 
aircraft  center  of  gravity  and  the  force  on  the  gear 
from  these  sources  is  already  accounted  for,  the  strut 
displacements  from  static  load  sources  are  neglected. 
Instead,  only  the  dynamic  strut  displacements  are  consi- 
dered. Thus  the  simulation  begins  with  the  aircraft  in 
a three  point  contact  with  the  runway,  with  a sink  rate, 
and  the  dynamic  strut  displacements  at  equilibrium. 


To  model  each  strut  it  is  important  to  take  into  con- 
sideration the  geometry  of  the  aircraft  as  seen  in  Fig  6. 


Fig.  6.  Strut  Compression  Mechanics 


The  compression  of  the  main  gear  strut,  for  example,  will 
result  from  two  separate  motions;  a vertical  displacement 
of  the  aircraft  center  of  gravity  will  change  the  strut 
displacement,  and  a pitching  rotation  about  the  center  of 
gravity  will  change  the  strut  displacement.  The  same 
applies  for  the  nose  gear  with  the  exception  of  a sign 
change  for  a rotational  displacement.  Expressing  this 
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in  equation  form,  one  arrives  at 


Zc  = X»0  + Z ft. 


zsng  ='XB0  + Z ft. 


where  ZSmg  and  ZSng  are  the  strut  displacements,  0 is  the 
angular  displacement  of  the  aircraft  from  horizontal  and 
Z is  the  vertical  displacement  of  the  aircraft  center  of 
gravity. 

In  a similar  fashion  one  arrives  at  the  rate  of  change 
of  displacement  for  each  strut. 

Z.  = X.0  + Z ft. /sec. 
smg  A 

ZSn^  = -X^0  + Z ft. /sec. 

Scaling  these  equations  for  the  computer  requires 


an  estimate  of 

the  maximum  value  of 

each  variable,  as 

summarized  in 

Table  I . 

Table  I 

■ 

Variable 

Max  Value 

Computer  Variable 

Z 

10  ft 

[Z/10] 

• 

Z 

100  ft/sec 

[Z/100] 

e 

0.25  rad 

[0/0.25] 

1 ® 

2.5  rad/sec 

[0/2.5] 

| 

1 
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Variable 

Max  Value 

Computer  Variable 

2 Sing 

1 ft 

IW1) 

• 

Zsmg 

10  ft/sec 

I2smg/10l 

2sng 

10  ft 

[ZSng/10] 

Zsng 

50  ft/sec 

fZsng/50] 

Arranging  these  equations  in  terms  of  computer 

variables 

leads  to 

^mg/1]  = XA[0/O.25)O.25  + [2/10)10 

ft. 

[ZSng/10]10  = -XB [0/0 .25)0.25  + [Z/10J10 

ft. 

[Zsmg/10110  = XA [ 0/ 2 . 5 J 2 . 5 + [Z/100J100 

ft. /sec. 

[ZSng/50]50  = -XB [0/2. 5)2.5  + [Z/100]100 

ft ./sec . 

Dividing  through  by  the  coefficients  of  the  left 
side,  substituting  for  XA  and  XBr  and  simplifying  leads  to 


[ZSmg/l]  = 0.8125  [0/0.25) 

+ 10[Z/10] 

ft. 

[ZSng/10J  = -0.500610/0.25) 

+ [Z/10J 

ft. 

[Zsmg/10]  = 0.8125 [0/2. 5)  + 

10[Z/100] 

ft. /sec 

tZSng/50]  = -1. OOl[0/2. 5]  + 

2 [Z/100] 

ft. /sec 

• • 

The  variables  0,  0,  Z and  Z are  determined  by  solving  the 

three  equations  of  motion  of  the  aircraft  and  the  static 

• • • . 

load  distribution  equation,  simultaneously,  for  X,  0,  and  Z. 


24 


I 


Once  the  respective  strut  displacements  and  rates  of  dis- 
placement are  known,  the  dynamic  force  contribution  to 
the  near  can  be  determined.  If  one  assumes  that  the  maxi- 
mum dynamic  main  gear  load  is  100,000  lbs.  and  the  maximum 
dynamic  nose  gear  load  is  50,000  lbs.,  one  arrives  at 

lFd*g/l°0.°00U°°.0°°  =,  Kdnig[ZSmg/I0]10  + KSmg[ZSmg/l] 

IFcl„g/50><100]50,000  - Kdng[ZSng/50150  + Ks„g  IzSng/101 10 

To  complete  the  above  equations  it  is  necessary  to 
determine  the  spring-rate  and  damping  coefficient.  The 
spring  rate  coefficients  may  be  determined  by  linearizing 
the  load  vs.  displacement  schedule  for  each  strut  about 
the  static  equilibrium  displacement  value.  For  the  F-4E 
the  main  gear  and  nose  gear  spring  rates  are  respectively 
120,000  lbs. /ft.  and  19,500  lbs. /ft.  (Ref  8:V2,5).  The 
viscous  damping  coefficients  for  each  main  gear  and  the 
nose  gear  as  determined  by  the  manufacturer  are  respectively 
5100  lb. -sec. /ft.  and  1170  lb. -sec. /ft.  (Ref  15:18).  Con- 
sidering that  there  are  two  main  gear  struts,  the  above 
equations  can  be  simplified  by  dividing  through  with  the 
left  side  coefficients. 

tpdm/1(>51  " 1-020  [Zsmg/10]  + 2.4|ZSmg/U 
|Fdng/5X104l  = 1.170 [Zsng/501  + 3.9(Zsng/10] 
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To  determine  the  static  force  on  the  landiny  gear  it 


L*. 


is  necessary  to  sum  the  external  vertical  forces  on  the 
aircraft  and  proportion  the  load  according  to  the  formula 


F 


smg 


Xa 

(xA  + xB) 


IF 


z 


Xd 


sng  (XA  + XB) 


IF. 


Tire/Wheel  Dynamics . With  the  completion  of  the 
aerodynamics,  engine  dynamics,  and  vertical  strut  dyna- 
mics, the  equations  of  motion  of  the  aircraft  are  readily 
solved,  if  one  assumes  a constant  value  for  the  coefficient  of 
friction  between  the  tire  and  runway.  Unfortunately  this 
coefficient  of  friction  is  dependent  on  many  factors  in- 
cluding runway  surface  condition,  wheel  speed  and  aircraft 

i 

speed.  For  a given  tire  and  runway  condition,  however,  a 
characteristic  nonlinear  relationship  between  the  degree 
of  wheel  slip  and  coefficient  of  friction  will  always  hold 
(Ref  17:1).  This  relationship  is  known  as  a p vs.  a curve, 
where  p is  the  coefficient  of  friction  and  a is  the  ratio 

of  (X  - Vw)  to  X.  Recall  that  X is  the  true  airspeed  or 

* 

ground  speed  in  this  no-wind  case  and  vw  is  the  wheel  speed 
or  Vw  + R w,  where  Vw  is  the  wheel  speed  component  due 
to  fore-aft  strut  oscillation,  where  Rr  is  the  tire  rolling 
radius,  and  u is  the  wheel  angular  velocity.  To  determine 
the  wheel  slip  ratio,  a,  it  is  necessary  to  find  the 
wheel  speed,  Vw,  since  X has  already  been  determined  for 
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the  aerodynamics  model.  If  the  fore-aft  strut  oscilla- 
tion speed,  Vws,  is  determined  by  the  horizontal  strut 
dynamics  model  and  the  tire  rolling  radius,  Rr , is  assumed 
to  be  a constant,  then  the  wheel  angular  velocity,  w,  is 
the  only  unknown.  To  determine  w it  is  helpful  to  recall 
the  formula 


E (torques)  = IW<I> 

where  Iw  is  the  mass  moment  of  inertia  about  the  rolling 
axis  of  the  wheel.  For  a constant  mass,  constant  geometry 
wheel,  Iw  is  a constant.  To  determine  u it  is  only  neces- 
sary to  know  the  torques  acting  on  the  wheel.  For  a braked 
aircraft  tire,  assuming  the  normal  forces  act  through  the 
center  of  the  wheel,  the  only  torques  are  those  caused  by 
the  ground  on  the  tire  and  the  brakes  on  the  wheel.  Summing 
moments  about  the  center  of  mass  of  the  wheel,  one  obtains 

Iwu  = (Ground  Torque)  - (Brake  Torque) 

If  the  brake  torque  is  determined  in  the  brake  dynamics 
model,  the  only  unknown  is  the  ground  torque.  Working  with 
a single  main  gear  wheel,  the  torque  can  be  represented  by 

(Ground  Torque)  » 1/2  (yFZjng)  Rr 

where  l/2(yFZing)  is  the  load  on  one  main  gear  strut/tire. 
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Note,  that  in  a round-about  way,  it  takes  u to  be  able  to 
calculate  y.  On  the  analog  computer,  an  initial  value  of 
wheel  speed  will  provide  an  initial  value  for  y,  therefore, 
the  loop  is  self  perpetuating  once  begun.  To  scale  this 
for  the  computer,  once  again,  requires  an  estimate  of  the 
maximum  values  for  the  variables  as  summarized  in  Table  II. 


Variable 


Table  II 
Maximum  Value 
10*  rad/sec^ 

10*  rad/sec 
105  lbs 
25,000  ft- lbs 
10*  ft/sec 


Computer  Variable 
[u/10*] 
tw/10*] 

[BT/25000 ] 
[Vw/10*] 


Substituting  for  the  constants,  the  equation  of  motion  for 
the  wheel  is 


(2 .71) [ ®/10*]  10*  = (0.5)  [yFZjng/105]105(1.16) 

-IBT/25000] 25000 


Dividing  through  by  the  coefficients  of  the  left  side  and 
simplifying,  one  arrives  at 


[ w,/10*l  = 2.140 [yF_  /105]  - 0.9225  [BT/25,000] 


By  integrating  w , mis  known.  With  u>,the  wheel  slip  ratio 
is  known.  To  obtain  a it  is  necessary  to  scale  (x  - Vw) . 
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Recalling  that  X had  been  originally  scaled  to  280  ft. /sec. 
to  be  compatible  with  (Vw/104],  it  may  now  be  rescaled  to 
a maximum  value  of  10,000  by  multiplying  by  0.0280. 
Therefore,  one  arrives  at 

[a/1]  = (0.0280 [X/10000]  - [Vw/10000])  * 0. 0280  [X/10000] 

For  the  particular  tire  on  the  F-4E  the  u vs.  o relationship 
is  plotted  at  several  velocities  and  several  runway  surface 
conditions  from  dry  to  flooded  (Ref  17) . The  general  form 
of  this  curve  is  shown  in  Fig.  7. 


Fig.  7.  Coefficient  of  Braking  Friction  vs.  Wheel  Slip  Ratio 

By  programming  a Variable  Diode  Function  Generator  (VDFG) 
on  the  analog  computer,  the  y vs.  a relationship  can  be  re- 
presented by  a ten  straight  line  segment  approximation.  To 
correct  this  curve  for  a velocity  other  than  the  one  it  is 
plotted  against,  one  may  assume  a linear  relationship  of  Ay  to 
AX  about  a central  velocity.  For  example,  the  central  velocity 
selected  for  this  problem  was  75  knots  or  approximately  mid-way 

between  the  maximum  touchdown  velocity  and  a stop.  The  VDFG  was 

• 

programmed  with  y vs.  o at  75  knots.  To  determine  the  Ay  vs.  AX 
relationship,  the  maximum  values  of  y at  0 knots  and  150  knots 


wore  compared  against  the  maximum  value  of  p at  75  knots. 
A median  value  of  Ap  over  the  two  75  knot  spans  was  se- 
lected and  also  programmed.  The  resulting  equation  was 


[Ap/1]  = (0.1849  - 8.684 [Vw/104])  x [p/1] 

Combining  the  p and  Ap  produces  a coefficient  of  friction 
as  a function  of  runway  surface  condition,  wheel  speed, 
and  aircraft  speed. 

Brake  Dynamics . The  brake  dynamics  model  must  be 
designed  to  develop  the  brake  torque  required  by  the  tire/ 
wheel  dynamics  model.  To  do  this  several  aspects  of  the 
brake  hydraulic  system  must  be  considered.  Assuming 
that  the  aircraft  brake  hydraulic  pump  delivers  a constant 
pressure  source  of  hydraulic  fluid,  brake  torque  is  deve- 
loped dependent  on  antiskid  system  inputs,  pressure  lag, 
the  pressure  to  torque  relationship,  brake  fade,  and 
torque  peaking.  The  antiskid  system  inputs  will  be  deve- 
loped in  the  antiskid  dynamics  model  and  are  assumed  to 
be  known  at  this  point.  The  pressure  lag  characteristic 
of  the  F44E  system  is  modeled  as  a first  order  lag  trans- 
fer function  with  a 0.1  second  time  constant  as  is  done 
in  McDonnell's  Motion  Base  Simulator  (Ref  3241) . The 
Laplace  transform  for  this  lag  can  be  written  as  10/(s+10) 
which  is  readily  programmed  on  the  analog  computer.  The 
pressure  to  torque  relationship  for  the  F-4E  brake  as 
plotted  in  Air  Force  Aeronautical  System  Division/ENFL 
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TM-73-1  takes  the  form  as  shown  in  Fig.  8 


I 


Brake 
Torque 

Brake  Pressure 

Fip.  8.  Brake  Pressure  Vs.  Brake  Torque 

By  programming  a Variable  Diode  Function  Generator  this  re- 
lationship between  brake  torque  and  brake  pressure  can  be 
stored  in  the  analog  computer.  To  model  the  brake  fade 
phenomena,  the  technique  used  in  Boeing's  Brake  Control 
Simulator  was  used  (Ref  15:129).  By  conservation  of  energy 

Eb  a Wb  = jBT  d©  * J BTo>  dt 

where  is  the  energy  absorbed  by  the  brake  heat  sink 
(assuming  no  energy  losses  to  convection,  radiation, 
conduction,  etc.),  Wb  is  the  work  done  by  the  brake,  0 is 
the  angular  displacement  of  the  wheel,  u>  is  the  wheel  angu- 
lar velocity,  and  t is  time.  The  conversion  from  energy 
to  brake  temperature  is  required  so  that  fade  initiation 
can  begin  once  a threshhold  temperature  is  reached.  This 
is  done  by  dividing  the  energy  absorbed  by  the  brake  by 
the  mass  of  the  brake  heat  sink  and  the  specific  heat  of 
the  heat  sink  material.  Thus,  the  equation  becomes 

31 


L 


Boeing  in  their  brake  control  simulator  suggested  a 
brake  fade  initiation  threshhold  of  1000°  F above  ambient. 
Using  this  value  as  a switch,  above  which  brake  torque 
fades,  then  a formula  for  ABT  as  a function  of  temperature 
would  describe  the  fade  phenomena.  For  this  model  a 
value  of  ABTp=  0 . 5BT/1000  F was  used. 

The  torque  peaking  effect  or  the  tendency  of  brakes 
to  increase  brake  torque  at  an  increasing  rate  as  a stop 
is  approached,  was  patterned  after  the  Boeing  Brake 
Control  Simulator  model.  Boeing's  analysis  of  brake 
torque  histories  showed  a pattern  as  depicted  in  Fig  9. 


Fig.  9.  Brake  Torque  Vs.  Time 

The  curve  shows  a quasi-steady-state  level  of  torque  with 


a sharp  increase  at  the  end  of  the  ground  roll.  The 
characteristic  value  of  this  torque  peak  at  the  end  of 
the  ground  roll  was  found  to  be,  on  the  average,  25%  of 
the  quasi-steady-state  torque  value.  Using  this  value  and 
a wheel  speed  of  90  ft. /sec.  as  a threshhold,  the  average 
initiation  value  for  F-4E  type  brakes,  a formula  for  torque 
peaking  is  arrived  at 

AP.Ttp  = (0.0463  - 5.140 [Vw/104] ) * [BT/25000] 


32 


By  combining  the  antiskid  command,  pressure  lag,  pressure 
to  torque  relationship,  brake  fade,  and  torque  peaking, 
the  brake  dynamics  model  is  complete. 

I 

Antiskid  Dynamics ■ The  Hytrol  Mark  II  antiskid  system 
used  on  the  F-4E,  basically,  measures  the  wheel  speed,  dif- 
ferentiates it  with  respect  to  time,  and  compares  this 
value  of  wheel  deceleration  against  reference  values  to 
determine  when  and  how  much  hydraulic  pressure  to  damp 
(Ref  14:55).  A deceleration  below  approximately  20  ft/sec2, 
will  result  in  full  system  pressure  being  delivered.  Above 
20  ft/sec2  and  below  140  ft/sec2,  a fraction  of  the  brake 
pressure  proportioned  to  (Vw  - 20)/ (140-20)  will  be  dumped. 

Vw  in  this  case  is  the  wheel  deceleration.  Above  1'40  ft/sec^ 
full  brake  pressure  is  dumped.  Modeling  this  on  the  analog 
computer,  however,  requires  a value  for  Vw  to  be  able  to 
calculate  Vw.  To  do  this  a feedback  loop  must  be  created. 

To  accomplish 


(Vw  - 20)  = Vw  ~ 20 

(140  -20)  120 


• 

Vw  must  be  scaled  and  a division  operation  performed.  The 
scaling  causes  the  rumerator  value  of  this  expression  to  be 
on  the  order  of  10“^,  a value  of  marginal  accuracy  on  the 
analog  computer.  When  division  is  performed,  the  inaccuracy 
is  multiplied  and  when  this  value  is  processed  through  the 
feedback  loop,  further  distortion  occurs.  Hence,  the  wheel 
deceleration  is  lost  in  the  background  noise  of  the  computer. 
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Instead  of  the  exact  reproduction  of  the  antiskid  system, 
a facsimile  was  tried  where  the  results  of  a real  system 
were  modeled.  McDonnell,  on  their  Motion  Base  Simulator, 
derived  a formula  for  the  on-off  action  of  the  pressure  dump 
valve  based  on  a square  wave.  The  on  portion  of  the  square 
wave  was  a fraction  of  the  total  period  of  the  wave  and  de- 
pendent on  the  wheel  speed.  The  total  period,  in  turn,  was 
a second  function  of  the  wheel  speed  (Ref  3:41). 

(On  Time)  = (1  - 0.001075  Vw)  (Period) 

(Period)  = (1.2  - 0.005  Vw) 

To  model  this  on  the  analog  computer  would  require  an  ex- 
panded memory  capability  or  a restartable  time  reference 
to  keep  track  of  the  elapsed  time  of  each  period  in  the 
square  wave.  Neither  of  these  were  feasible  without  re- 
sorting to  a digital-analog  hybrid  system.  In  an  attempt 
to  keep  the  system  model  totally  analog,  a third  procedure 
was  attempted  where,  once  again,  the  results  of  the  real 
system  were  modeled.  In  this  approach  the  wheel  speed  was 
measured  against  the  aircraft  speed,  and  controlled  to  re- 
main in  the  region  X > Vw  > (X  - 60).  The  60  ft. /sec. 
value  was  arrived  at  by  observing  traces  of  wheel  speed 
vs.  time  for  real  aircraft  landings,  wherein  the  wheel 
speed  cycled  between  the  aircraft  velocity  and  this  value. 
When  the  aircraft  reached  a speed  25  ft. /sec.  the  antiskid 
system  was  by-passed  and  full  brake  pressure  was  delivered; 
otherwise  a 5 ton  would  have  been  aporoached  asymptotically. 
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To  keep  the  wheel  speed  within  the  physical  range  of 
0 < Vw  < X,  comparator  circuits  had  to  be  installed  that  set 
Vw  * 0 or  Vw  = X depending  on  the  actual  calculated  value  of  u>. 
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Fore-Aft  Strut  Dynamics . The  phenomena  that  caused  a 
tire  to  develop  a coefficient  of  braking  friction  is  the 
occurence  of  a velocity  differential  between  the  tire  and 
contacting  surface.  Since  the  relative  velocity  of  the  tire 
sensed  by  the  runway  is  at  issue,  all  components  of  wheel 
velocity  must  be  considered.  By  limiting  the  landing  air- 
craft problem  to  three  degrees  of  freedom,  the  lateral  or 
side  to  side  motion  of  the  tire  is  eliminated,  thus,  leaving 
only  rotational  and  translational  components.  The  rotational 
component  of  wheel  speed  was  addressed  by  the  tire/wheel 
dynamics  model;  the  translational  motion  remains. 

If  all  the  structures  between  the  wheel  and  aircraft 
center  of  gravity  were  rigid,  the  wheel  translational  velo- 
city would  be  identical  to  the  ground  velocity  and  no  speed 
differential  would  exist.  Since  the  landing  gear  is  not  rigid, 

i 

its  equation  of  motion  must  be  considered  to  determine  this  component. 

By  taking  the  landing  gear  strut  to  be  a cantilever 
beam,  it  is  known  that  any  load  in  the  form  of  a force 
or  a torque  will  result  in  a linear  and  an  angular  dis- 
placement due  to  the  elastic  effects.  However,  in  the 
case  of  a cantilever  beam  where  the  deflections  are  small, 
such  as  a landing  gear  strut,  a reasonable  assumption  can 
be  made  to  uncouple  the  linear  displacement  from  the 
angular.  For  the  F-4E,  the  fore-aft  spring  cons- 
tant is  30-1, COO  lbs. /ft.  at  static  equilibrium  strut 
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extension.  If  the  maximum  strut  load  is  50,000  lbs.  and 
the  maximum  tire  coefficient  of  friction  is  0.7,  the  max- 
imum braking  force  is  35,000  lbs.  Such  a force  would 
cause  a deflection  of  1.38  inches.  When  considered  over 
the  length  of  the  strut,  3.91  feet,  this  results  in  an 
angular  deflection  of  1.68  degrees,  a negligible  amount. 
Therefore,  a mass  spring  damper  analogy  can  be  made  where 
the  mass  of  the  unsprung  strut  and  wheel  is  concentrated 
in  the  wheel  and  the  strut  serves  as  a spring-damper  com- 
bination for  linear  movement  in  the  fore-aft  direction  as 
depicted  in  Fig*  10. 
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Fig.  10.  Fore- Aft  Strut  Analogy 


If  the  damping  coefficient  for  the  F-4E  is  200  lb. -sec./ 
ft. , the  equation  of  motion  becomes 

MgXg  - -200Xg  -304 , 800XS  -0.5yFZmg 

Substituting  scaled  variables  and  dividing  through  by  the 
coefficients  of  the  left,  the  equation  takes  the  form 


[Xg/50000]  = -0.6275 [Xg/2000]  -9. 562 [Xg/20] 


-0.0784  [yFZing/105] 
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Since  a positive  strut  velocity  exacerbates  the  velocity 
differential  between  tire  and  runway,  the  Xs  is  actually 
subtracted  from  the  wheel  velocity  to  produce  the  relative 
velocity  of  the  tire  to  runway. 

Total  System  Integration.  The  completion  of  the  sub- 
system models  now  allows  the  three  equations  of  motion  and 
the  static  load  distribution  equation  to  be  solved  in  the 
form  of  the  aircraft  dynamics  model.  Starting  with  the 
summation  of  forces  in  the  x-direction  as  depicted  in 
Fig.  4,  one  obtains 

•• 

ZFX  = MX  = Tcos*t  - D - yFZmg  - WrFZng  - Ddc 

where  M is  the  mass  of  the  aircraft,  T is  the  thrust,  #T 

is  the  thrust  inclination  angle,  D is  the  aircraft  drag, 

y is  the  braking  coefficient  of  friction,  FZjng  is  the 

total  main  gear  load,  yr  is  the  rolling  coefficient  of 

friction,  FZng  is  the  total  nose  gear  load,  and  Ddc  is  the 

drag  of  the  drag  chute. 

•• 

If  X is  scaled  to  have  a maximum  value  of  40  ft. /sec./ 
sec.  and  yr  is  equal  to  a constant,  0.025,  then  the  scaled 
equation  takes  the  form 


M[X/40]  40  = [T/l] (0.9958)  - [D/1]  - [yFz  /105] 

mg 

- 0.025[F,ng/105)  - [Ddc/1] 


Substituting  for  [T/l]  , [D/1],  and  [D<jc/1]  with  the  values 
developed  previously  and  dividing  by  the  coefficients  of 
the  left  side,  the  equation  takes  the  form 

fX/40]  = 23,4014-  - 2-  ‘ 7A.  [X/280 ] - 144.5106  fx/28012 

M M M 

- ilHlyF.  /105]  - 62  -5l°-°-[Fz  /105] 

M 2 mg  M n9 

- 231.5875 f»/2ani2 

M 

If  X has  a maximum  value  of  8000  ft.,  the  equations  for 
X and  X are 

[X/280]  * 0.1429  MX/140]  dt  + X(0) 

[X/8000]  * 0.0350  / [X/280]  dt  + X(0) 


The  equation  representing  the  summation  of  forces  in 
the  z-direction  is  actually  two  separate  equations  due  to 
the  distinction  between  quasi-static  and  dynamic  gear 
loads.  For  the  quasi-static  gear  loads,  the  equation  can 
be  stated  as 


Fsmg  + Fsng 


mg 


L - Tsin4>T 


where  Fsmg  and  FSng  are  the  main  and  nose  gear  static 
loads,  respectively,  and  mg  is  the  weight  of  the  aircraft. 
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Scaling  FSmg  + FSng  to  a max  value  of  100,000  lbs.  and 
substituting  scaled  equations  for  L and  T,  one  arrives  at 

I(F_  + F_  )/105]  = 3.2174  x 10'4(M)  - 0 . 1344 [X/280] 2 
amg  ^ng 

- 0.0009  + 0.00010  [X/280] 

By  specifying  the  center  of  gravity  location  the  proportion 
of  the  total  gear  load  going  to  the  nose  gear  is  determined 
by  the  ratio  XB  f (XA  + XB)  where  XB  is  the  horizontal  dis- 
tance from  the  center  of  gravity  to  the  nose  gear  and 
XA  is  the  horizontal  distance  from  the  main  gear  to  the 
center  of  gravity.  For  a 0.33c  center  of  gravity  this 
ratio  is  equal  to  0.1396  implying  that  the  ratio  of  main 
gear  load  to  total  static  load  is  0.8604.  Therefore  the 
static  load  distribution  equation  becomes 

t(F=mg  + Fs„g>/105)  - 0-1396t(Fsmg  ♦ FSng)/10b 

+ 0.8604  l(FSmg  ♦ F3ng)/105) 
tFSmg/1<)51  = 0.8604tFsmg  + FSng)/10s1 
tFs„c/50.000]  = O.U96[(FSmg  + FSng)/105] 

For  the  dynamic  gear  loads  the  equation  of  motion  can 
be  stated  as 
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• 1 '* 
.ill 


*luCJ 


+ F, 


ng 


If  Z has  a riaxir.um  value  of  200  ft. /sec. the  scaled 

values  for  Fj  and  are  substituted,  and  the  equa- 

amg  ang 

tion  is  divided  by  the  coefficients  of  the  left  side,  the 
equation  takes  the  form 


IZ/200]  = -^ptFdmg/105]  -i|0[Fdng/50,000] 

Note  that  the  negative  signs  allow  a positive  force  to 
compress  the  strut. 

If  Z and  Z have  maximum  values  of  100  ft. /sec.  and 
10  ft.,  respectively,  then  equations  for  [Z/100]  and 
tZ/10]  can  be  written  as 


tZ/100]  = 2 JtZ/200]  dt  + Z(0) 
tZ/10]  = 10  jtZ/100]  dt  + Z (0) 

Summing  moments  about  the  center  of  gravity  of  the 
aircraft  as  depicted  in  Fig  4,  the  pitching  moment  equa- 
tion becomes 


**Mcg  = Iyy0  ~ Mp  + Ddczdc  + FZngXB  ~ wrFzngzng 
~ FzmgxA  " uFZmgZmg 


where  Iyy  is  the  aircraft  mass  moment  of  inertia  about 
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the  y-axis,  0 is  the  angular  acceleration,  flp  is  tlie 
pitching  moment  caused  by  the  untrimmed  flight  controls, 

Zjjc  is  the  vertical  distance  from  the  drag  chute  attach- 
ment point  to  the  aircraft  center  of  gravity,  FZng  and 
FZjng  are  the  summation  of  static  and  dynamic  gear  loads 
on  the  nose  gear  and  main  gear,  respectively,  and  Zng  and 
zmg  are  identical  constants  representing  the  vertical 
distance  from  the  runway  to  the  aircraft  center  of  gravity. 

By  substituting  scaled  equivalents  for  terms  deve- 
loped in  the  subsystem  models,  by  scaling  0 to  a maximum 
value  of  25  rad. /sec. /sec. , and  by  substituting  values 
for  the  moment  arms,  the  equation  can  be  expressed  as 

Iyy[0/25]25  = 16641 . 691 [X/280] 2 + 31403 . 26 [X/280] 2 

+ 1, 001, 250[FZng/50, 000] 

- 7887.5 [FZng/50, 000]  - 325 , 000 [FZmg/105] 

- 631,000 tuFZmg/105] 

Dividing  through  by  the  coefficients  of  the  left  and  com- 
bining the  two  nose  gear  terms,  the  equation  takes  on  its 
final  form 

[0/25]  = i£iliil[X/280]2  + 12A6’130[X/280]2 

Iyy  lyy 

+ 39734’5[Fz ng/50,000]  - i!20°[Fz  /105] 

Jyy  Jyy 
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If  0 and  0 have  maximum  values  of  2.5  rad. /sec.  and 
0.25  rad.,  respectively,  then  the  equations  for  [0/2.5] 
and  [0/0.25]  are 

[0/2.5]  = 10  I [0/25]  dt  + 0(0) 

[0/0.25 ]=  10  | [0/2 .5]  dt  + 0(0) 

The  analog  comouter  wiring  diagram  for  the  total  system 
integration  to  this  point  is  depicted  in  Fig.  11.  The 
variables  are  all  scaled  to  their  respective  maximum 
values.  Triangles  represent  summers,  and  triangles  with 
rectangles  on  the  left  side  represent  integrators.  The 
sign  of  a variable  is  reversed  when  passing  through  the 
summers  and  integrators  to  account  for  apparent  reversed 
signs.  The  oblong  circles  represent  attenuators  through 
which  the  variables  are  scaled  and  proportioned  according 
to  their  appropriate  coefficients.  The  computer  program 
in  this  figure  would  be  sufficient  to  calculate  main  gear 
loads  if  a constant  coefficient  of  friction  between  the 
tire  and  runway  were  assumed. 

By  referring  to  wiring  flow  charts,  the  process  of 
producing  a coefficient  of  friction  based  on  wheel  speed, 
aircraft  velocity,  and  runway  surface  condition  is  more 
easily  grasped.  In  Fig.  12,  page  44,  the  top  row  of 
summers  and  integrator  produces  the  wheel  speed  given  the 
ground  torque  and  brake  torque.  The  second  row,  consisting 
of  comparators  and  relays,  keeps  the  wheel  speed  within 
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: 


the  range  0 < Vw  < X so  that  the  wheel  slip  velocity  is 
only  in  one  direction.  The  reason  for  this  becomes 
evident  in  the  third  row  where  the  wheel  slip  velocity  is 
divided  by  the  sum  of  the  aircraft  velocity  and  the  fore- 
aft  strut  component  to  the  wheel  velocity  to  produce  o, 
the  wheel  slip  ratio.  The  sign  of  o must  always  be  posi- 
tive to  conform  to  the  physical  realities.  Continuing 
in  the  third  row,  the  Variable  Diode  Function  Generator 
produces  p'  based  on  the  median  velocity  of  75  knots. 

In  the  fourth  row,  p is  modified  for  any  wheel  velocity 
other  than  75  knots.  The  values  of  the  coefficients  in 
the  fourth  row  correspond  to  dry  runway  conditions. 

In  Fig.  13,  the  brake  dynamics  model  is 
depicted.  The  top  row,  consisting  of  a Variable  Diode 
Function  Generator,  provides  the  brake  torque  component 
produced  by  brake  pressure.  The  second  row,  relying  on 
a summer,  a comparator,  a relay,  and  a multiplication 
circuit,  produces  the  brake  torque  component  resulting  from 
operating  in  the  torque  peaking  velocity  range.  The 
third  row,  consisting  of  a multiplication  circuit,  an 
integrator,  a summer,  and  a comparator/relay  combination, 
produces  the  brake  torque  component  due  to  brake  fading. 

The  sum  of  the  three  brake  torque  components  results  in 
the  brake  torque. 

In  Fig  14  page  47,  the  antiskid  and  fore-aft  strut 
dynamics  models  are  depicted.  The  antiskid  model  uses  a 
series  of  two  comparators  and  relays  to  test  the  wheel 
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slip  speed  for  compliance  with  the  range  0 < (X  - Vw)  < 60 

ft. /sec.  If  so,  full  brake  pressure  is  delivered  to  the 

brakes.  If  not,  zero  brake  pressure  is  delivered  to  the 

• 

brakes.  A test  is  then  made  on  the  aircraft  velocity,  X, 
using  the  third  comparator  and  relay  to  check  if  the  air- 
craft  velocity  is  in  the  range  0 <_  X _<  25  ft. /sec.  If 
so,  full  brake  pressure  is  delivered  to  the  brakes.  In 
the  second  row,  the  touchdown  protection  feature  of  the 
antiskid  system  is  modeled.  By  comparing  the  time  elapsed 
against  a two-second  reference  value,  a two-second  delay 
on  brake  application  is  provided.  Furthermore,  this 
delay  allows  the  wheel  to  spin-up  to  a near-aircraft  velo- 
city before  braking  commences.  In  the  third  row,  the 
brake  pressure,  BP',  is  modified  by  the  10  * (S  + 10) 
transfer  function  to  provide  a first  order  lag.  In  the 

fourth  row,  the  fore-aft  strut  dynamics  are  modeled  to 
• • 

provide  Xs,  where  Xs  equals  VWg,  the  wheel  speed  component 
due  to  strut  oscillation.  Physical  constants  pertinent 
to  the  F-4E  are  listed  in  Appendix  A.  Source  data  neces- 
sary to  complete  the  model  is  listed  in  Appendix  B. 

Data  Acquisition 

To  test  the  aircraft  dynamics  model,  initial  condi- 
tions  were  placed  on  the  -X  integrator  corresponding  to 
the  touchdown  velocity,  on  the  -Z  integrator  correspond- 
ing to  the  sink  rate  and  on  the  to  integrator  corresponding 
to  a value  that  would  cause  the  wheel  speed  to  be  near 
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the  aircraft  velocity  two  seconds  after  touchdown.  All 
other  integrators  had  initial  conditions  of  zero.  In- 
herent in  this  procedure  is  the  assumption  that  the  wheel 
speed  will,  indeed,  spin-up  in  the  allotted  time.  For 
all  but  flooded  and  icy  runways  this  is  a good  assumption 
(Ref  14:218).  It  is  important  to  note  that  this  proce- 
dure is,  not  only,  a necessary  boundary  condition  on  the 
mathematics  for  the  model  to  simulate  a real  landing, 
but  it  is  a requirement  for  the  actual  Mark  II  antiskid 
system  to  ensure  proper  operation  (Ref  14:7).  For  a 
given  gross  weight  the  touchdown  speed  is  determined  by 
the  flight  manual  (Ref  13:B18).  The  sink  rate  is  arbi- 
trarily determined  to  provide  the  maximum  main  gear 
strut  deflection  without  entering  into  a bounce  condition, 
where  the  main  gear  would  lose  contact  with  the  runway. 

Test  runs  were  made  for  wet  and  dry  runway  y vs.  a 
data  (Ref  17).  For  each  runway  condition,  test  runs  were 
made  at  gross  weights  of  30,000  lbs.,  35,000  lbs., 

40,000  lbs.,  and  45,000  lbs.  At  each  gross  weight,  the 
mass  moment  of  inertia  of  the  aircraft  was  adjusted  accord- 
ingly. Actual  values  for  the  test  run  variables  are 
listed  in  Table  III* 

Data  Analysis 

To  determine  the  validity  of  the  model,  the  gear 
loadings  should  conform  to  actual  test  results  of  gear 
loadings.  In  the  absence  of  these,  however,  it  is 


Table  III 


Independent  Variables  for  Each  Data  Run 


Iyy  (ft-lb-SeC^) 

30,000 

1.074X105 

*o 

(knots) 

127 

Zo 

(ft/sec) 

5.40 

Qo 

(rad/sec) 

200 

M 

(slugs) 

932.43 

Gross  Weight  (lbs) 


35,000 

40,000 

45,000 

1.086X105 

1.098X105 

1.110X10 

137 

147 

155 

5.50 

5.80 

5.94 

270 

370 

380 

1087.83 

1243.24 

1398.64 

necessary  to  validate  the  test  results  circumstantially. 

Assuming  that  the  basis  of  the  model  is  sufficient 
to  cover  all  the  phenomena  occuring  in  a landing,  it  was 
shown  that,  primarily,  the  aerodynamics,  engine  dynamics 
and  vertical  strut  dynamics  produce  the  gear  load,  while 
the  tire/wheel  dynamics,  brake  dynamics,  antiskid  dyna- 
mics, and  fore-aft  strut  dynamics  are  required  to  deter- 
mine the  stopping  distance.  If  the  coefficient  of  fric- 
tion, u,  is  accurate,  the  gear  loads  should  be  accurate 
within  the  scope  of  the  assumptions  since  the  aerodynamics, 
engine  dynamics , and  landing  gear  spring  constants  and 
damping  coefficients  are  well  documented. 

To  determine  if  the  coefficient  of  friction,  y,  is 
accurate  consider  the  latter  four  subsystems.  In  the 
case  of  the  tire/wheel  dynamics  model,  the  actual  p vs.  o 
data  is  well  documented.  To  determine  if  a,  the  wheel 
slip  ratio,  is  accurate  one  must  know  if  the  aircraft 


velocity,  fore-aft  strut  velocity,  and  wheel  velocity  are 
accurate.  Aircraft  velocity  will  be  accurate  if  y,  once 
again,  is  accurate,  because  the  aerodynamics,  engine 
dynamics,  and  strut  dynamics  can  be  modeled  with  a high 
degree  of  certainty.  Since  fore-aft  strut  spring  cons- 
tants and  damping  coefficients  are  well  documented,  the 
fore-aft  strut  velocity  should  be  accurate.  Wheel  speed 
accuracy  can  then  be  traced  to  braking  action  accuracy. 
Since  the  brake  torque  vs.  brake  pressure  relationship 
is  well  documented,  the  accuracy  of  y can  be  traced  back 
to  the  accuracy  of  the  antiskid  system.  To  determine  if 
the  model  approximates  real  aircraft  antiskid  operation, 
traces  of  wheel  speed  vs.  time  and  stopping  distance  will 
be  enlightening.  If  the  simulation  wheel  speed  wave 
form  approaches  the  actual  wheel  speed  wave  form  and  the 
stopping  distances  are  accurate,  then  the  model  should 


be  accurate 


III.  Results  and  Conclusions 


\ 

t- 

\ 


System  checkout  revealed  that  fore-aft  strut  dyna- 
mics contributions  to  the  wheel  velocity  occured  only  at 
the  instant  of  touchdown  and  were  damped  out  within  the 
first  two  seconds  thereafter.  Considering  the  fact  that 
the  brake  system  is  locked-out  by  the  touchdown  protec- 
tion circuit  until  two  seconds  after  touchdown,  the 
effect  of  the  fore-aft  strut  model  was  deemed  negligible 
and  bypassed.  Also,  system  checkout  revealed  that  at 
no  time  did  the  brake  fade  phenomena  occur.  Considering 
that  the  brake  energy  absorbed  is  characteristically  30% 
to  50%  of  that  for  a rejected  take-off,  for  which  the 
brakes  are  designed,  this  result  is  quite  reasonable. 

Comparison  of  the  plots  for  both  wet  and  dry  runway 
runs,  contained  in  Appendix  B,  shows  a marked  difference 
in  the  effect  of  a more  slippery  surface  on  the  initial 
deceleration  at  touchdown  and  on  the  stopping  distances. 
For  the  eight  combinations  of  gross  weight  and  runway 
surface  condition,  the  stopping  distances  are  assembled 
in  Table  IVand  compared  to  flight  test  stopping  distance 
as  extracted  from  the  flight  manual  (Ref  13:B8**2)  . 

In  general,  the  model  brakes  to  a stop  in  a shorter 
distance  than  the  aircraft.  The  maximum  stopping  dis- 
tance differences  are  6.25%  on  dry  runways  for  the  40,000 
lbs.  aircraft  and  6.10%  on  wet  runways  for  the  45,000 
' lbs.  aircraft. 
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Table  IV 

Stopping  Distance  of  Model  and  Aircraft  (ft.) 


r 


Model  Aircraft 


Gross  Weight 

I 

Dr* 

Wet 

Dry 

Wet 

30,000 

2300 

4590 

2400 

4600 

35,000 

2640 

5480 

2800 

5600 

40,000 

3000 

6270 

3200 

6500 

45,000 

3390 

6760 

3500 

7200 

I Figure  17,  page  64,  shows  the  decelerating  effect  of 

brakes  and  drag  chute  quite  prominently  at  two  and  four 
seconds,  respectively.  When  zero  aircraft  velocity 
occurs,  instabilities  result  from  the  tire/wheel  dynamics 
model  attempting  to  divide  by  zero  and  cause  erratic 
plotting  pen  movement.  The  plots  of  landing  gear  load 
and  angular  motion  also  demonstrate  a lurching  tendency 
as  the  stop  is  reached. 

With  the  possible  exception  of  the  antiskid  system,  the 
subsystem  models  behave  well.  The  aerodynamics,  as  indica- 
ted indirectly  in  plots  of  aircraft  velocity  and  stopping 
distance,  show  the  characteristic  parabolic  curve  as  the  ve- 
locity decreases.  Furthermore,  the  assumption  of  constant 
angle  of  attack  on  the  runway  appears  well  justified  after 
considering  the  brief  influence  of  the  aircraft  pitch  angle 
in  the  0 vs.  time  plots.  The  engine  dynamics,  though  not 
plotted,  follow  the  predicted  values.  The  vertical  strut 
loads  correspond  to  both  the  dynamic  strut  loads  and 
the  quasi-static  strut  loads  as  depicted  in  the  main 
and  nose  gear  load  plots.  In  these  it  is  important  to 
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note  the  damping  of  the  dynamic  loads  and  the  gentle 
increase  in  static  loads  as  the  lift  decays.  Plots  of  o, 
li,  V , and  X when  compared  show  that  y is  accurate  for  a 
given  a,  and  a is  accurate  for  the  given  aircraft  velocity 
and  wheel  velocity.  The  acceleration  of  the  wheel, 
based  on  ground  torque  and  brake  torque  is  responsive  to 
changes  in  the  coefficient  of  friction,  y,  thus  indicating 
that  the  tire/wheel  dynamics  model  produces  reasonable 
results.  Observation  of  the  brake  pressure  and  brake 
torque  plots  indicate  strong  cyclic  dependence  on  the 
antiskid  system  commands.  A strong  decrease  in  brake 
pressure  output  as  the  wheel  velocity  approaches  zero 
velocity  illustrates  the  high  frequency  of  the  antiskid 
response  and  its  ability  to  control  wheel  velocity  in 
this  region. 

It  is  important  to  note  that  due  to  scaling  require- 
ments and  the  gain  factor  required  to  plot  wheel  velo- 
city, the  noise  level  becomes  a factor  in  tending  to  mask 
the  wheel  velocity  wave  form.  Figure  15  indicates  the 
nature  of  the  real  aircraft  wheel  velocity  wave  form  for 
dry  and  wet  runways. 

For  the  wet  runway  model,  the  pronounced  "scalloping" 
effect  is  not  realized,  although  both  the  real  aircraft 
wheel  velocity  and  the  model  wheel  velocity  tend  to  ride 
50-60  ft. /sec.  less  than  the  aircraft  velocity.  The  lack 
of  the  "scalloping"  effect  indicates  that  the  model  reacts 
much  quicker  than  the  real  aircraft  to  wheel  velocity 


***0  = 256  ft./soc. 
Ref.  15:67,  68 


Time  (sec) 


Wet  Runway- 


Time  (sec) 


Fig.  15.  Aircraft  Wheel  Speed  vs.  Time 


deviations.  This  is  to  be  expected  when  recalling  that 
the  first  two  attempts  at  modeling  the  antiskid  system 
ran  into  difficulty  when  confronted  with  an  inability  of 
the  analog  computer  to  model  a variable  wave  form  with 
a specified  minimum  response  time.  Since  the  dry  runway 
wheel  velocity  plots  are  comparable,  a conclusion  might 
be  drawn  to  search  for  a better  wet  runway  antiskid  sys- 
tem, although  stopping  distances  are  generally  better 
correlated  than  dry  runway  distances.  An  ideal  substi- 
tute would  .be  an  actual  antiskid  system  in  the  loop, 
although  all  that  is  really  needed  is  an  ability  to 
command  a specified  "off"  time  during  which  full  brake 
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pressure  is  dumped.  An  antiskid  system  of  this  nature 
would  result  in  a more  realistic  wet  runway  wheel  velo- 
city plot  and  tend  to  lengthen  the  ground  roll  slightly. 

Altogether  the  model  gives  reasonable  stopping  per- 
formance and  system  simulation.  Assuming  that  the  simu- 
lated stopping  distances  are  within  tolerance  and  the 
antiskid  model  is  reasonably  accurate,  a schedule  of  load 
and  velocity  vs.  distance  may  be  constructed  from  cross 
plots  of  main  gear  load,  velocity,  and  stopping  distance 
vs.  time.  This  alternative  schedule  is  shown  in  Fig  16. 


Runway  Di stance  ( ft . ) 


Velocity 


100  3000 
Runway  Distance  (ft.) 


Fir.  16.  Derived  Main  Gear  Load/Velocity 


Schedule 


IV.  Recommendations 


All  indications  are  that  an  improved  antiskid  system 
will  increase  system  accuracy  considerably.  The  most 
likely  method  of  simulation  to  achieve  this  is  a hybrid 
analog-digital  computer  system.  Such  a system  would 
eliminate  the  weakness  of  the  analog  system  by  providing 
precise  control  over  the  brake  pressure-off  time  and 
increase  the  resemblance  of  the  two  wet  runway  wheel 
velocity  vs.  time  plots. 

Another  recommendation  is  to  place  an  actual  braked 
tire/wheel  with  antiskid  in  the  loop.  By  operating  such 
a system  on  a dynamometer  wheel  wherein  the  exact  load  on 
the  tire  is  controlled  by  the  analog  program,  highly 
accurate  performance  should  result.  Furthermore,  such  a 
procedure  would  open  up  increased  analog  capacity  to 
model  the  noh-linear,  dual-chamber  nature  of  the 
F-4E  main  gear  strut. 

Finally,  adoption  of  this  system  simulation  procedure 
for  tire  and  landing  gear  development  is  urged.  Total 
analog  simulation  of  the  aircraft  dynamics,  despite  draw- 
backs to  the  antiskid  model,  gives  a reasonably  accurate 
schedule  of  gear  loading  vs.  velocity  which,  when  compared 
to  the  current  military  specification,  is  considerably 
more  realistic.  Further  improvements  in  the  antiskid  model 
by  resorting  to  digital-analog  models  or  by  using  actual 
antiskid  systems  in  the  circuitry,  should  provide  highly 
realistic  results. 
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Appendix  A 


Physical  Constants 


Contained  herein  are  the  specific  values  of  terms 
that  were  assumed  to  be  constant.  Included  are  the 
physical  constants  and  geometric  measurements  of  the  F-4E 
that  are  essential  to  solving  the  equations  of  motion. 


Physical  Constants 


Roman  Svmhols 


Greek  Symbols 


c 

CD 

cdbasic 

Cd<Jc 

cL 

CLBASIC 

Cm 

CmBASIC 

cPb 

g 

Kdmg 
Kdmgh 
Kdng 
Ksmg 
^Smgh 


KS 

Mb 

Ms 

Rr 

S 


ng 


XA 

Xb 

Zdc 

zmg 

zng 


16.04  ft 

0.117 

0.140 

0.1875 

0.272 

0.195 

0.021 

0.025 

0.11  BTU/lbm°F 
32.174  ft/sec2  * 
2.71  ft-lb-sec2 
5100  lb-sec/ft 
200  lb-sec/ft 
1170  lb-sec/ft 
120,000  lb/ft 
304,800  lb/ft 
19,500  lb/ft 
3.85  slugs 
12.75  slugs 
1.16  ft 
530  ft2 

3.25  ft 

20.025  ft 
3.39  ft 
6.31  ft 
6.31  ft 


a 

&CDge 

^CLge 

ACMge 

Vi 

V 

p 


-1° 

-0.023 

0.077 

-0.004 

0.025 

3.14159 

0.002378  lb-sec/ft2 
5.25° 
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Appendix  B 


Modeled  Aircraft  Performance  vs . Time 

Contained  herein  are  the  test  results  of  the  system 
simulation  in  5000  lb.  increments  of  aircraft  gross  weight 
for  dry  and  wet  runway  conditions.  For  each  gross  weight 
and  runway  condition  there  are  six  figures  containing 
plots  of  the  following  variables  versus  time. 

First  figure:  X;  X;  X 

«•  • 

Second  figure:  0;  0;  0 

Third  figure:  Z;  Z;  Z 

Fourth  figure:  BP';  BP;  BT 

Fifth  figure:  o;  v;  Vw 

Sixth  figure?  FZmg;  FZng;  l»FZmg 


ESS 


Fig.  18.  0,  0,  and  0 vs.  Time  for  30,000  lb.  F-UE,  Dry  Runway 


(ft. /sec.) 


19.  Z,  Z,  and  Z vs  Time  for  30,000  lb.  F-hS,  Dry  Runway 


SES 


■m 


-30 

-20 

-10 


Fig.  25.  Z,  Z,  and  Z vs.  Time  for  35,000  lb.  F-UE,  Dry  Runway 


Fig.  28.  F , F , andAF  vs.  Time  for  35,000  lb.  F-UE 
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Fig.  39. O’ ,/a,  and  V vs.  Time  for  h^tnno  lb.  F-hE,  Dry  Runway 


•0.025 


BP,  and  BT  vs.  Time  for  30,000  lb.  F-LE,  Wet  Runway 


Pig.  UU.  CT  ,m  , and  vs.  Time  for  30,000  lb.  F-lt5,  Vfet  Runway 


Pig.  U6.  F , F , and/«F_  vs.  Time  for  30,000  lb.  F-US,  Wet  Runway 


(ft. /see 


Fig.  h9 . Z,  Z,  and  Z vs.  Time  for  35,000  lb.  F-UE,  Vet  Runway 


Fig.  SO.  BP',  BP,  and  BT.vs.  Time  for  35,000  lb.  F-UE,  Wet  Runway 


Fig.  5l.<T  ,Mt  and  V vs.  Tine  for  35,000  lb.  F-hH,  Wet  Runwa; 


Fig.  5U.  0,  0,  and  0 vs.  Time  for  L0,000  lb.  F-US,  Wet  Runw 


Wet  Runway 


lb.  F-UE,  Wet  Runway 


and  0 vs.  Time  for  U5,000  lb.  F-ljS,  Wet  Runway 
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Pig.  62.  BP',  BP,  and  BT  bs.  Time  for  LS,000  lb. 
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